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We study the spectral curves of dinier models on periodic Fisher graphs, obtained 
from a ferromagnetic Ising model on Z^. The spectral curve is defined by the zero locus 
of the determinant of a modified weighted adjacency matrix. We prove that either they 
are disjoint from the unit torus (T^ = {{z,w) : \z\ = 1, \w\ = 1}) or they intersect 
at a single real point. 

1 Introduction 

u 

(-H In this paper we study the spectral curve of periodic, 2-dimensional Fisher graphs, either 

finite in on direction, and periodic in the other direction (cylindrical graph); or bi-periodic 
^ and obtained from a bi-periodic, ferromagnetic Ising model on (toroidal graph). To an 

I— —I edge-weighed, cylindrical (resp. toroidal) Fisher graph, one associates its spectral curve 

P{z) = 0{resp.P{z,w) = 0). The real polynomial P{z){resp.P{z,w)) defining the spectral 
'j>! curve arises as the characteristic polynomial of the Kasteleyn operator in the dimer model. 

*0 The study of spectral curve for periodic Fisher graphs (which is non-bipartite), is 

inspired by the work of Kenyon, Okounkov and Sheffield [ KO061 IKOS06j . They prove 
fT) that the spectral curve of bipartite dimer models with positive edge weights is always a 

QQ real curve of a special type, namely it is a Harnack curve. This implies many qualitative 

and quantitative results about the behavior of bipartite dimer models related to the phase 
transition. 

• • The Fisher graph we consider in this paper is a graph with each vertex of the honey- 

• ^ comb lattice replaced by a triangle, see Figure 1. 

^ A planar graph is one which can be embedded into the plane such that edges can 

intersect only at vertices. Fix an embedding of a planar graph G. A clockwise-odd 
orientation of G is an orientation of the edges such that for each face (except the outer 
face) an odd number of edges pointing along it when traversed clockwise. For a planar 
graph, such an orientation always exists |Kas67j . The Kasteleyn matrix corresponding to 
such a graph is a |F(G)| x |1/(G)| skew-symmetric matrix K defined by 

V and n — )■ V 

Ku,v = { —W{uv) if M ~ t> and u <— v 

else. 

where W{uv) > is the weight associated to the edge uv. 

Now let G be a Z^-periodic planar graph. By this we mean G is embedded in the plane 
so that translations in Z^ act by weight-preserving isomorphisms of G which map each 
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Figure 1: Fisher graph 

edge to an edge with the same weight. Let Gn be the quotient graph G/{n'L x nZ). It is 
a finite graph on a torus. Let 7a;^„(resp. 7j,,n) be a path in the dual graph of G„ winding 
once around the torus horizontally(resp. vertically). Let Eh (resp. Ey) be the set of edges 
crossed by 72,- (resp. 7^^). We give a crossing oriGntcitioii for the toroidal graph G^i as 
follows. We orient all the edges of G„ except for those in Eh U Ey. This is possible since 
no other edges are crossing. Then we orient the edges of Eh as if Ey did not exist. Again 
this is possible since G — Ey is planar. To complete the orientation, we also orient the 
edges of Ey as if Eh did not exist. 

Let Ki be a Kasteleyn matrix for the graph Gi. Given any parameters z,w^ we con- 
struct a matrix K{z., w) as follows. Let "jx^i and 'jy^i be the paths introduced as above. 
Multiply Ku^y by z if the orientation on that edge is from u to v, and multiply K^^v by 
^ if the orientation is from v to u, and similarly for w on "/y. Define the characteris- 
tic polynomial P{z,w) = det K{z,w). The spectral curve is defined to be the locus 
P{z,w) = 0. 

We also discuss cylindrical graphs in this paper. Assume we have a Fisher graph, 
periodic in the z direction, finite in the w direction, by removing all i/;-edges and ^ edges. 
An example of such a graph is shown in Figure 2. 

Assume it is m x n, meaning that it has width m with respect to z, but period n 
with respect to w. We embed this graph into a cylinder. Let K(z) be the corresponding 
weighted adjacency matrix with an orientation shown in Figure 2. K{z) is obtained from 
K{z,'w) by giving all edges crossed by 7^ weight 0. Let P{z) = det K{z). 

Kenyon and Okounkov ([KO06]) proved that the spectral curves of any periodic bipar- 
tite graph with positive edge weights are Harnack curves \Mik\ IMROlj , whose intersection 
with Tx,y = {{z,w) : \z\ = e^, \w\ = e^} can only be: i. no intersection; ii. a pair of conju- 
gate points, each of which is of multiplicity 1; Hi. a single real zero of multiplicity 2 (real 
node). Since the Fisher graph we consider in this paper is non-bipartite, no previous result 
is known. We prove the following theorems 

Theorem 1.1. Consider a positively-weighted cylindrical Fisher graph, which is periodic 
in one direction, and finite on the other direction. Its quotient graph under the translation 
can be embedded into a cylinder, as illustrated in Figure 2. The intersection of P{z) = 
and T is either empty or a single real point. 



2 



Figure 2: Fisher graph on a cy Under 



In Figure 1, each edge of the graph is either a side of the triangle (triangular edge), 
or an edge connecting different triangles (non-triangle edge). We call all the horizontal 
non-triangle edges in Figure 1, o-edges. 

Theorem 1.2. Consider a bi-periodic Fisher Graph, given edge-orientations as illustrated 
in Figure 1. Assume all the triangular edges and a-edges have weights 1, and all the other 
edges have weights in (0,1), then 

1. the spectral curve P{z,w) = is a Harnack curve; 

2. the only possible intersection of P{z,w) = with the unit torus = {[z,w) : \z\ = 
l,\w\ = 1} is a single real point of multiplicity 2. 

Our result is very promising. Firstly, it leads to important properties of the dimer model 
on cylindrical graphs, such as weak convergence of Boltzmann measures and convergence 
rate of correlations for an infinite, periodic graph with finite width along one direction, see 
Proposition 3.4. Secondly, since the dimer model on the Fisher graph are closely related 
to the Ising model |MW| and the vertex model [Lilj . our result leads to a quantitative 
characterization of the critical temperature of the arbitrary periodic ferromagnetic, two- 
dimensional Ising model, as the solution of an algebraic equation. |Li2) . 




The outline of the paper is as follows. In Section 2, we explain the connection between 
the Kasteleyn operator, characteristic polynomial, spectral curve with the dimer model. 
In section 3, prove Theorem 1.1, as well as a phase transition characterized by the decay 
rate of edge-edge correlation for the dimer model on the cylindrical Fisher graph, resulting 
from Theorem 1.1. In section 4, we prove Theorem 1.2. The proof consists of 2 critical 
components, one is that P{z, vu) > for any {z, vo) € T^, the other is an explicit correspon- 
dence between the Kasteleyn operator on the Fisher graph, and the Kasteleyn operator in 
the square-octagon lattice, introduced in |Dubl2] . 
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2 Background 



A perfect matching, or a dimer cover, of a graph is a collection of edges with the 
property that each vertex is incident to exactly one edge. A graph is bipartite if the 
vertices can be 2-colored, that is, colored black and white so that black vertices are adjacent 
only to white vertices and vice versa. 

To a weighted finite graph G = {V, E, W), the weight W : E ^ is a function from 
the set of edges to positive real numbers. We define a probability measure, called the 
Boltzmann measure /x with sample space the set of dimer covers. Namely, for a dimer 
cover D 

^^{D) = ^ n ^(^) 

eeD 

where the product is over all edges present in D, and Z is a normalizing constant called 
the partition function, defined to be 

D ee-D 

the sum over all dimer configurations of G. 

If we change the weight function W by multiplying the edge weights of all edges incident 
to a single vertex v by the same constant, the probability measure defined above does not 
change. So we define two weight functions W, W' to be gauge equivalent if one can be 
obtained from the other by a sequence of such multiplications. 

The key objects used to obtain explicit expressions for the dimer model are Kasteleyn 
matrices. They are weighted, oriented adjacency matrices of the graph G defined as in 
Page 1. 

It is known |Kas6H IKas671 ITesOOl IKQSOGj that for a planar graph with a clockwise odd 
orientation, the partition function of dimers satisfies 

Z = Vdet K. 

Given a Fisher graph with an orientation as illustrated in Figure 1, the quotient graph 
can be embedded into an n x n torus. When n is even, if we reverse the orientations of 
all the edges crossed by 7^. and all the edges crossed by 7^, the resulting orientation is a 
crossing orientation. For 9,t G {0, 1}, given the orientation as in Figure 1, let Kn^ be the 
Kasteleyn matrix Kn in which the weights of edges in Eu are multiplied by (—1)^, and 
those in Ey are multiplied by (—1)^. Using the result proved in |TesOO] . we can derive that 
when n is even, the partition function Z„ of the graph Gn is 

Zn = l\- Pf(i^r) + Pf(^r) + Pf(i^r) + Pf(^^')|- 

Let Em = {ei = uivi, Cm = UmVm} be a subset of edges of G„. Kenyon |Ken97j 
proved that the probability of these edges occurring in a dimer configuration of Gn with 
respect to the Boltzmann measure P„ is 

Pnieu em) = ^^='Yz^^'"'\ - mK')E^^ + + Pf{K')E^ + Fml')E^J 

where E^ = V{Gn) \ {ui, vi, Um, Vm}, and {K^'^)e^ is the submatrix of Kfj' whose lines 
and columns are indexed by E!^. 
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The asymptotic behavior of Z„ when n is large is an interesting subject. One important 
concept is the partition function per fundamental domain, which is defined to be 

lim ^logZ„. 

n— >oo n 

Gauge equivalent dimer weights give the same spectral curve. That is because after 
Gauge transformation, the determinant is scaled by a nonzero constant, thus not changing 
the locus of P{z,w). 

A formula for enlarging the fundamental domain is proved in [CKPOll IKOSOGj . Let 
Pn{z,w) be the characteristic polynomial of Gn, and Pi{z,w) be the characteristic poly- 
nomial of Gi, then 

Pn{z,w) = n n 



3 Graph on a Cylinder 
3.1 Spectral Curve 

Proof of Theorem 1.1 Without loss of generality, assume the period n, or the circum- 
ference of the cylinder, is even. Assume P{z) = has a non-real zeros zq G T^. Assume 

^o = e^"o-, ao G (0,1)U(1,2) 

We classify all the real numbers in (0, 1) U (1, 2) into 3 types 

1. ao = I where p, q are positive integers with no common factors, p is odd 

2. ao is irrational 

3. oo = I where p, q are positive integers with no common factors, p is even 
First let us consider Case 1 and Case 2. There exists a sequence ik G N, such that 

lim Zq'' = — 1 

fc— >oo 

In other words, if we assume Zq*" = e^'^"*'^ where a^ G [0, 2), then 

lim Ofc = 1 

k—^oo 

According to the formula of enlarging the fundamental domain, 

Pizo) = 

implies 

Pi,{z'on) = o 

Since the cylindrical graph is actually planar, if we reverse the orientation of all the edges 
crossed by 'jx, we get a clockwise-odd orientation, given that n is even. Hence 
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where Zi^, is the partition function of dimer configurations of the cyhnder with circumfer- 
ence nl}; and height m. Therefore we have 



1 = hm 



^4(-l) 



(1) 



= lim 



fc— >-oo Z"^ 




(2) 



'4 




[7r(afc 



(2t)! 



(3) 



Since 



^4(-)= E ^^'''^(-. + ^) 



0<j<m 



where '^'^ is the signed sum of loop configurations winding exactly j times around the 
cylinder (see Lemma 2.1), we have 



To estimate * |, let us divide the m x ^^n cylinder into ^^n, m x 1 layers. Connecting 
edges between two layers may be occupied once, unoccupied, or occupied twice(appear as 
doubled edges). Choose one layer Lq, we construct an equivalent class of loop configura- 
tions. Two loop configurations are equivalent if they differ from each other only in Lq, and 
coincide on all the other layers and boundary edges of Lq. 

We claim that for any given equivalent class, there is at least one configuration including 
only even loops. To see that we choose an arbitrary configuration in that equivalent class 
including odd loops. Let T be the set of all triangles in Lq belonging to some loop crossing 
Lq. We choose an odd loop s\. Choose an triangle Ai G si H Lq. Moving along Lq until 
we find a triangle A2 belonging to a different non-planar odd loop S2- This is always 
possible because by Lemma 2.2, any loop configuration including odd loops always has 
an even number of odd loops, and all the odd loops are winding once along the cylinder. 
Change path through Ai. Starting from Ai, we change the doubled edge configuration to 
alternating edges along a path in Lq, and change the path through all triangles, belonging 
to even loops, or si, between Ai and A2, then we change the path through A2. For any 
even loop between A2 and A2, we change paths through an even number of triangles of 
that even loop, hence it is still an even loop. However, for s\ and S2, we change paths 
through an odd number of triangles, then both of them become even. We can continue this 
process until we eliminate all the odd loops, because there are always an even number of 




l<j<m 



and 
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odd loops in the configuration. An example of such an path change process is illustrated 
in the following figures. 

Since we have at most 2™ different configurations in each equivalent class, and each 
equivalent class has at least one configuration including only even loops, we have 



# of loop configurations including odd loops 
# of even loop configurations 



< 2^ 



We have a finite number of different edge weights, each of which is positive, hence the 
quotient of any two weights is bounded by a constatn Ci, then 

|pj^fc)| < Partition of configurations including nonplanar odd loops 

l<j<m 

+Partition of configurations with only even loops 
< (2'"Cf™ + 1) Partition of even loop configurations < C^Zf 



As a result 



lim — 2" 



E 

t=l 



[7r(a,-l)]2*02*Pfc(e 



< lim — K- 
k^oo Zf ^ 
4- t=i 



{2t)\ de^i 
[rmriak - 1)]^* 



^2 



Since m is a constant, and limfc_>.oo ctfc = 1, we have 



lim > 

k^oo ^ — ' 



t=l 



[mnjak - 1)] 

mi 



2t 



-Co" 



which is a contradiction to ([3|). Hence for all ao in Case 1 and Case 2, P(e*"'"^) 7^ 0. 

Now let us consider oq in Case 3. From the argument above we derive that as long as 

|a - 1| < 5, 

where (5 is a small positive number depending only on m, 

P(e^"^) / 0. 

Consider ao = ^, where p, q have no common factors and p is even. We claim that as long 
as the denominator 1 is sufficiently large, P{zq) 7^ 0. To see that, if 

1 6 

there exists an integer k, such that 



-Arg[e''"^°''] - 1 



TT 



< s. 



Hence after enlarging the fundamental domain to an m x kn cylinder, we have 

Pkie'^'^'^n + 0, 
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hence 

Now we consider | > |, only finitely may q's satisfy this condition. Let ^ be a prime 
number satisfying £ > [|] + 1, then after enlarging the fundamental domain to an m X in 
cylinder, the corresponding 6 will not change because it depends only on m. For any 

. 1 3 
q 2 



we have 



Let £ai{ mod 2) = ao- Namely, 



6 1_ 
2^ £q' 



P + 2sq 

ai = — for s = 1,2,..., £ 

£q 

in reduced form. By the previous argument P(e"i*'^) 7^ 0, hence P^(e""'^^^^) 7^ 0. Hence 
when ao is in Case 3, Pe{zo) 7^ after enlarging the fundamental domain to m x n£, where 
i depends only on m, and is independent of ao- If P{zo) = 0, zq is not real, then we enlarge 
the fundamental domain to m x ni, where £ is a big prime number depending only on m, 
we derive that -P^(^o) — however this is impossible since Pe{z) = can have only real 
root on T. □ 

Corollary 3.1. For any non-real z £T, all eigenvalues of K{z) are of the form iXj,j = 
1, ...,6mn, where i is the imaginary unit. 3mn of the Xj 's are positive, the other Xj 's are 
negative. 

Proof From the definition of K{z), iK{z) is a Hermitian matrix. We claim that iK{z) 
has 3mn positive and Smn negative eigenvalues for any non-real z E T. In fact, for 
a planar graph with no z vertices, Kq is a anti-symmetric real matrix with eigenvalues 
zbiAj, 1 < J < 3mn — m, and deti^o > 0, because \PfKo\ is the partition function of 
dimer configurations with positive edge weights, and deti^o = {PfK^f' By the previous 
theorem, every time we add a pair of boundary vertices connected by a 2; edge, we have a 
anti-Hermitian matrix K^+i of order Qmn — 2m-\-2{r + 1) with as a principal minor. By 
induction hypothesis iK^. has the same number of positive and negative eigenvalues, the 
interlacing theorem implies that iK^+i has at least 2>mn — m + r positive and 3mn — m + r 
negative eigenvalues. By previous theorem, det K^+i > for non-real z, so the other two 
eigenvalues of iKr+i can only be one positive and one negative. □ 



Lemma 3.2. If P{1) = 0, then = 0. Let z = e^^ , then ^^^''^'^ 

2=1 



dz 



0. For 



generic choice of edge weights, —q^ 



z=l 



Proof. We prove the result for derivatives with respect to z, the derivatives with respect 
to Q are very similar. Since 



p(.) = ;^p,(.^ + i: 



9 



we have 



dP 
dz 



2 = 1 



Z=l 



3 



2 = 1 



Each monomial in P(l) corresponds to a loop configuration including only even loops. 
Consider a configuration including 2k even loops, each of which winding exactly once 
around the cylinder. Except those non-contractible even loops, the rest of the graph is 
covered by doubled-edge configuration. Consider P(l) as a polynomial of edge weights. 
The monomial corresponding to that configuration has a coefficient 

Q I c\2k 

>^0,k = ±^ 

because each single loop can have two different orientations, corresponding to two terms in 
the expansion of the determinant. However, the coefficient of the monomial corresponding 
to the same configuration in ^ ^Jf^ ^ is 

S2,k = {2kf + 2k{2k - 2)2 + 2fc y2k-4f + ---+(^liy2k- 4kf 



Obviously, is a number depending on k. Therefor P(l) and ^^-^ 
each other. Let 



cannot divide 



2 = 1 



Wo 



eJeeE 



d'^P{z) 



0} 



2=1 



Wo = {{We)eeE : P{1) = 0} 



The intersection of W2 and Wq forms a proper subvariety of W. Hence if P{1) = 0, and 



we choose the edge weights generically, ^ q^^^ 



2=1 



□ 



Remark. Lemma 3.3 implies that if -P(l) = 0, for generic choice of edge weights, the 
intersection is of multiplicity 2. 



Proposition 3.3. Let 



1 



27r 



F{r) = / log P{re"')dd 
27r Jo 



If P{z) = does not intersect the unit torus T, F{r) is differentiable at r = 1; If P{z) = 
has a real zero of multiplicity 2 at z = 1, 



lim 



dF{r) 



r-->i+ dr r->i- dr 
Proof. Since P{z) is a Laurent polynomial in z, according to the Jensen's formula 



dF{r) 



iog|p(o)i^-|:.og(^)+i/ 



27r 
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where ai, • • • , a„ are the zeros of P in the interior of the disk {z : \z\ < r}. If P{z) has no 
zeros on the circle {z : |z| = r}, then 

dF ^ Fir + Ar) - Fir) 

-TT- = hm 

or Ar-!>o Ar 



Um -— log f ) 

ir^o Ar ^ V r 



k=l 
n 

r 

The proposition follows from substituting r = 1, and the fact that the intersection of 
P{z) = with T can only be a single real point. □ 

Remark. Jensen's formula implies that the curve P{z) = is Harnack if and only if 

dF(r) 

ro( lim - lim ) ^ < 2, 

r'^ro+ r— !>ro- Or 

for any rg > 0. In fact, if the height of the cylinder m < 3, we can always derive that the 
corresponding spectral curve is Harnack. To see that, first of all, being Harnack is a closed 
condition, for generic choice of edge weights, the intersection of P{z) = with \z\ = 1 
is at most two points(counting multiplities) . Without loss of generality, assume for some 
r(0 < r < 1), P{z) = intersects \z\ = r at 3 different points zi,Z2,Z3, then ^1,^2,^3, 
lie also on the intersection of P{z) = and \z\ = r, given that P{z) is a real-coeffient 
polynomial. Hence the intersection of P{z) = with \z\ = r is at least 4 points. Moreover, 
by symmetry, the reciprocal of those points are also roots of P{z) = 0, then P{z) = has 
at least 8 roots, which is a contradiction to the fact m < 3. 

3.2 Limit Measure of Cylindrical Approximation 

From the proof we know that all terms in Pmxni^ — 1) are positive, if n is even. We can 
always enlarging the fundamental domain in z direction without changing edge weights to 
get 

Pmx2n(-1) = Pmxnii)Pmxni-i) = ^™xnW = \P fmx2nK 

Therefore Pmxn{i) is the partition function of dimer configurations of the m x 2n cylinder 
graph. According to the formula of enlarging the fundamental domain, 



mx2l 



n{ 1) — Pmxn{z), 



,2i- 



we have 



lim T7logPmx2Zn(-l) = 1™ Yl ^OgPmxniz) (4) 



A;^4/ ° '""^'"^ ' A;^47r2/ 



1 f dz 

— I \ogPmxn{z)— (5) 



The convergence of the Riemann sums to the integral follows from the fact that the only 
possible zeros of P{z) on T is a single real node. ^ is defined to be the partition function 
per fundamental domain. 
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Any probability measure on the infinite banded graph, with depth m on one direction, 
and period n on the other direction, is determined by the probabihty of cyhndrical sets. 
Namely, we choose a finite number of edges ei, 62, • • • , arbitrarily, and the probabilities 

Pr(ei&e2& • • • fecfc) 

that 61,62, •• • ,6^ occur in the dimer configuration simultaneously for all finite edge sets 
determines the probability measure. We consider the measures on the infinite graph as 
weak limits of measures on cylindrical graphs. First of all, we prove a lemma about the 
entries of the inverse Kasteleyn matrix using the cylindrical approximation. 



Lemma 3.4. 



/— >-oo 



27r 



p.v. / Z' 



k,-k^ coJoclorK,ny,„{s,,s,,){z) dz 



IZ 



Proof. To that end, we construct a transition matrix S to make S ^K^-^2inS block diag- 
onal, with each block corresponding to a m x n quotient graph. Define 



C _ (-pi 6mri 1 6mn 1 6mn \ 



where 



■ i7r(2j + l)(2fe+l) 

e 4i 



s = i 
s^t 



then 



S-^K, 



inx2ln'- 



V 







Since S ^ = ^S*, we have 

^mx2lni~^){kv,Sv),{kw,Sw) 



21-1 







Kmxn{s 



i(4i-l)7i 



= -y 

21 ^ 



CofactorKmxn{Sv, Sw) i(2j+l)T^(kv-kyj) 



j=0 



(2j + l)i-n 



where is the index of the fundamental domain for vertex v and Sy is the index of vertex 



in the fundamental domain. 

If Pmxn{z) has no zero on T, we have 



^|™^-^mx2«n( ^){kv,Sv),{kw,Sw) ~ 2tt jj 



1 f ky-kyj(^ofactorKmxn{sv,Sw){z) dz 



Pmxn{z) 



IZ 



(6) 



If 1 is an order-2 zero of Pmxn{z), let 

cofactorKmxn{s 



Qiz) = z' 



Pmxn{z) 
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Since deti^(l) is an anti-symmetric real matrix of even order, and non-invertible, the 
dimension of its null space is non-zero and even. Hence AdjK{l) is a zero matrix, and 1 is 
at least a zero of order 1 for cofactorKmxnis^, Sn]){z). Then in a neighborhood of 1, 



, , Res^=iQ{z) 
Q{z) = + R{z), 



z-1 



where R{z) is analytic at 1. 



1 



lim iiT^io, ( — „ wi. o ^ = lim lim — 



Z— >oo 



+ 



E 



21- 



For the second term. 



lim lim — 



E (2j+l)i^ 
Q{e 2' ) 



lim — 
(5^0+ 27r 



TT 2 

27r-(5 



27r 



dz 



iz 



For the first and the third term 



lim lim — 

(5-5-0+ /-s>oo 21 



1 

lim lim — 

S^0+ l^cxj 21 



\ReSz=iQ{z) + ^ i?es|2|<i^^ 



(8) 

(9) 
(10) 

(11) 



\0<i< 

/ 



E + E 



Since i?(2:) is analytic in a neighborhood of 1, 



(2j + l)»7I 

Q(e 2i 



ReSz=iQ{z) 

(2j + l)i7r 

, e 2i — 1 



(2j + l)»7r , 
21 



lim lim — 

5^0+ ;^oo 21 



( 



\ 



E + E 



Vo<i<f-2 



(2j + l)i7r . 

-R(e 2! 



Moreover, 



lim lim —ReSz=\Q{z) 



E + E 



2Z- 



l-^<i<2ny 



(2j + l)t7i 

e 2i 



lim lim —ReSz=iQ(z) 

5^0+Z^oo 2/ ^ 7 ^ 



1 



+ 



- lim —ReSz=iQ{z) 
(5^0+ 21 



15 1 
7 + 2 



(2j + l)i7r ' (2j + l)i7i 

1 V e 21 — 1 e 2i 



And the theorem follows from ([6]), 0, ([TT]) , ([12]) , ([Ts]). 



(12) 



(13) 

(14) 

(15) 

□ 
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Proposition 3.5. Using a large cylinder to approximate the infinite periodic banded graph, 
we derive that the weak limit of probability measures of the dimer model exists. The prob- 
ability of a cylindrical set under this limit measure is 



k 



-1 Ui Vi ■ ■ ■ Uk Vk 



Pr(ei&e2&- • -fecfc) = lim T] We^\ det K .,, 

l^oo-^^ M/ mxzin y ... 

j=l V 

where We^ is the weight of ej, and Uj,Vj are the two vertices of ej 
Proof. On a finite cylindrical graph of m x 2ln 

Z, 



Pr(ei&e2&---&efc) 



z 



where Z is the partition function of dimer configurations on that graph, and ^ei,--- .e^. is the 
partition function of dimer configurations for which ei, • • • Cfc appear simultaneously. Since 

Z^ = detKm>c2ln{-l) (16) 
k 

-^ei,-,efc = Y{wejC0fact0rK^-^2ln,iuuVi,-,Uk,Vk)(.-'^) (17) 



(17) follows from the fact that if originally we have a clockwise-odd orientation, we still 
have a clockwise-odd orientation when removing edges and ending vertices, while keeping 
the orientation on the rest of the graph. The proposition follows from Jacobi's formula for 
the determinant of minor matrices. □ 

We consider two edges ei and 62 with weight xi and X2 on m x 2ln cylinder, and 
compute the covariance. 

-P?'mx2Zn(ei&e2) - Prmx2ln{ei)Prmx2ln{e2) (18) 



-X1X2 Jdeti^„i,2Z„ ( 2^ (-1) • detK^i,2/n [ll Zl) 



(20) 



= XiX2{\Kj^2lni'"l,Wi)Kj^2ln('"2,W2)+Kj^2lnivi^^2)Kj^^^^^ (21) 
-^mx2in(^'l'^'2)i^„x2te(^l'«^2)| " l^„x 2Zn("l ' ^1 )^mx 2«n(«2 , ^x;2) | ) (22) 

In order to compute the covariance as / oo, we only need to compute the entries of 

Then we have the following proposition 

Proposition 3.6. Consider the dimer model on a Fisher graph, embedded into an mx In 
cylinder, as illustrated in Figure 4- Let the circumference of the cylinder go to infinity, 
i.e. / — )• oo, and keep the height of the cylinder unchanged. Consider two edges ei and 
62- As |ei — 62! 00; if P{z) = does not intersect T, the edge-edge correlation decays 
exponentially ; if P{z) = has a node at 1, the edge-edge correlation tends to a constant. 



Proof. The theorem follows from formula ( 22 ) , and the estimates of the entries of inverse 



Kasteleyn matrix. By (11), the second term goes to as |ei — 62] — )■ 00, the first term is 0, 
if no zero exists on T, the first term is a nonzero constant if the spectral curve has a real 
node at 1. □ 
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Example 3.7. (1x2 cylindrical graph) Assume we have a dimer model on a Fisher graph 
embedded into an infinite cylinder of height 1 and period 2. One period of the graph is 
illustrated in Fig 5. 




Figure 5: 1 x 2 cylindrical graph 
The characteristic polynomial of the model is 

P{z) = (61622; — aia2)(^^ — 0102). 

The probability that an ai-edge occurs is 

01^2 I f ^ c^-^ f ^ dz \ 

Prlaij = — — [p.v. / — ^^^+p.v. / 7-r 

47r Jj aia2 - ^/^z h 0102 - ^ \f^z J 

{1 if 0102 > 6162 
if 0102 < 6162 
\ if 0102 = 6162 

P{z) = has a real node at 1 if and only if aia2 = 6162. At the critical case, the covariance 
of an ai edge and a 62 edge, as their distance goes to infinity, is 

Pr{ai&ib2) - Pr{ai)Pr{b2) = ^-\ = -\ 

4 Graph on a Torus 

4.1 Combinatorial and Analytic Properties 

To compute the characteristic polynomial P(z, w) of the Fisher graph, we give an orienta- 
tion to edges as illustrated in Figure 1. 

Lemma 4.1. The characteristic polynomial P(z,w) for a periodic Fisher graph with period 
{m, 0) and (0, n) is a Laurent polynomial of the following form: 

p(.,^) = ^p,,(.V + J-) 

hj 

where {i,j) are integral points of the Newton polygon with vertices (±m, 0), (0, ±n), (±m, ^n): 

Proof. Let p = 6mn. By definition, 

P{z, w) = det K{z, w)= X • • • X kp,i^. 

ll 
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(— m, n) 




(-m,0) 



(m,0) ^ 



(m, — n) 



Figure 6: Newton polygon 



Here a is the number of even cycles of the permutation (ii,...,ip). The sum is over all 
possible permutations of p elements. Each term of P{z, w) corresponds to an oriented loop 
configuration occupying each vertex exactly twice. For the graph Gi with m 2-cdges and 
n it;-edges, P{z,w) is a Laurent polynomial with leading terms z"^ ,w"' . z^w^ 
corresponds to loops of homology class P{z,w) is symmetric with respect to z'^w^ 

and . That is because for each term of z^w^ , if we reverse the orientation of all loops, 
we get a term of with coefficients of the same absolute value, corresponding to the 
product of weights of edges included in the configuration. The sign of the term is multiplied 
by (-l)f = 1. 

To show that all the powers lie in the polygon, we multiply all the 6-edges by z 
(or ^), and all the c-edges by w{or ^), according to their orientation. This way the corre- 
sponding characteristic polynomial becomes P{z"', w"^). Let (i, j) be a power of monomial 
in P(z^,w^). At each triangle, all the possible contributions of local configurations to 
the power of the monomial can only be (0, 0), (0, ±1), (±1, 0), (±1, ^1). Examples are 
illustrated in the following Figure 4. The left graph has two doubled edges, and the contri- 
bution to the power of the monomial is (0, 0), the right graph has a loop winding from the 
z-edge to the w-edge, and the contribution to the power of the monomial is (1, —1). 




Figure 7: local configurations 

Consider all the 2mn triangles. For each edge, we considered its contribution twice, 
since we counted it from both triangles it connected. Hence we have 

—2mn < 2i < 2mn 
—2nin < 2j < 2mn 
—2mn < 2(i + j) < 2mn 

Since i = ni, j = mj, and the Newton polygon N(P) is defined to be 

N{P) = convex hull{(i,j) G Z'^\z^w^ is a monomial in P(z,w)}, 

the lemma follows □ 

Lemma 4.2. For configurations with odd loops corresponding to a non-vanishing term in 
P{z, w), all odd loops have non-trivial homology, and the total number of odd loops is even. 
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Proof. For any loop configuration on the Fisher graph embedded on a torus, the number 
of odd loops is always even. That is because the total number of vertices are even, while 
odd loops always involve odd number of vertices. Any term in Pij including odd loops can 
appear only when odd loops have non-trivial homology. That is because for a contractible 
odd loop, we can reverse the orientation of that loop to negate the sign of that term. The 
term with reversed orientation on the odd loop cancels with the original term, because the 
homology class (0, 0) of the configurations is not changed given the odd loop is contractible. 

□ 

4.2 Generalized Fisher Correspondence 

Consider the Fisher graph obtained by replacing each vertex of the honeycomb lattice by a 
triangle. Assume all the triangle edges have weight 1, and all the non-triangle edges have 
positive weights not equal to 1. Furthermore, we assume that at each triangle, there is an 
even number of adjacent edges with weight less than 1. We introduce a generalized Fisher 
correspondence between the Ising model on the triangular lattice and the dimer model on 
the Fisher lattice as follows: 

1. If two adjacent spins have the same sign, and the dual edge of the Fisher lattice has 
weight strictly greater than 1, then the dual edge separating the two spins is present 
in the dimer configuration. 

2. If two adjacent spins have the same sign, and the dual edge of the Fisher lattice has 
weight strictly less than 1 , then the dual edge separating the two spins is not present 
in the dimer configuration. 

3. If two adjacent spins have the opposite sign, and the dual edge of the Fisher lattice 
has weight strictly greater than 1, then the dual edge separating the two spins is not 
present in the dimer configuration. 

4. If two adjacent spins have the opposite sign, and the dual edge of the Fisher lattice 
has weight strictly less than 1, then the dual edge separating the two spins is present 
in the dimer configuration. If at each triangle, an even number of incident edges have 
weight less than 1, we change the configuration on an even number of incident edges. 
As a result, the number of present edges incident to each triangle is still odd, which 
is a dimer configuration. This correspondence is 2-to-l since negating the spins at 
all vertices corresponds to the same dimer configuration. 

Around each triangle of the triangular lattice, we always have an even number of sign 
changes. If all the non-triangular edges have weight strictly greater than 1, then the number 
of present edges incident to each triangle of the Fisher graph is equal to the number of 
edges separating the same spins, which is odd. Whenever we have an edge with weight 
strictly less than 1, we change the configuration of that edge, according to the principle 
described above. Figure 9 is an example of the generalized Fisher correspondence given 
a > 1, & > 1, c < 1, d > 1, e < 1. 

Choose the interaction Jg associated to a bond as follows: 

Je = ^1 logiuel 

where Wg is the weight of the dual edge. J > corresponds to the ferromagnetic interaction. 
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Figure 8: Generalized Fisher Correspondence 
4.3 Duality Transformation 

We call a geometric figure built with a certain number of bonds of a lattice a closed 
polygon if at every lattice point only an even number of bonds occurs. It is clear that 
every configuration of "+" and "— " spins on a lattice can be associated a closed polygon 
of the dual lattice in the following way. A dual bond belongs to the polygon if it separates 
different spins and does not belong to the polygon if it separates equal spins. The same 
closed polygon is associated to two symmetric configurations in which the "+" and "— " 
spins are interchanged. 
Let 

Tmn be the quotient graph of the triangular lattice on the plane, as defined on Page 
4. Let Hmn be the dual graph of Tm„, Hmn is a honeycomb lattice which can be embedded 
into an m X n torus. Without loss of generality, assume both m and n are even. 

Define an Ising model on with interactions {^e}ee£;(T„„)- Assume the Ising model 
on Tjnn has partition function Zt^^j- Then Zt^^j can be written as, 

^T^n,/ = 2 n ®^P('^e) Yl Yl exp(-2Je) := 2 Y[ exp{-Je)ZF^^^Doo, 
e&E{Tmu) C*€S*oe€C* eeS(T^„) 

where Sqq is the set of closed polygon configurations of Hmn, with an even number of 
occupied bonds crossed by both and jy. The sum is over all configurations in ^qo- 
Similarly, we can define Sqi* {S^q, S^^) to be the set of closed polygon configurations of Hmn, 
with an even(odd,odd) number of occupied bonds crossed by 7a; , and an odd(even,odd) 
number of occupied bonds crossed by jy. 

Fmn is the Fisher graph obtained by replacing each vertex of Hmn by a triangle, with 
weights on all the non-triangle edges given by 



Let Zp^^^B be the partition function of dimer configurations on Fmn, a Fisher graph 
embedded into an m x n torus, with weights e^'^<= on edges of Hmn, and weight 1 on all the 
other edges. Then 

where 

ZF^„,De,^ = n ^M'^Je) E n exp(-2Je). 

For example, -^Fni^jDoi is the dimer partition function on Fmn with an even number of 
occupied edges crossed by ^x, and an odd number of occupied edges crossed by 7^. It also 
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corresponds to an Ising model which has the same configuration on the two boundaries 
parallel to 7y, and the opposite configurations on the two boundaries parallel to ^x- Similar 
results hold for all the Zf^^.Dq^^ 9,t G {0, 1}. 

On the other hand, if we consider the high temperature expansion of the Ising model 
on T^nj we have 



^Tmn,! = XI n exp(Jea„a^) 

o- e=uveE{T,n„) 

= X n (cosh Je + cr„(T^sinh Je) 
o- e=uveE{Tmn) 

= JJ cosh Je JJ (1 + (T„cr^ tanh Je) 

e=uveE{Tmn) °" e=UV<^E{Tmn) 

JJ cosh Je X n 

e=uveE{Tmn) CeS eeC 

where S is the set of all closed polygon configurations of Tmn- Let F„j„ be a Fisher graph 
embedded into an m x n torus, with weights tanh Jg on edges of G„, with each vertex of 
the triangular lattice replaced by a gadget, as illustrated in the following Figure. 





Figure 9: Fisher Correspondence 

There is an l-to-2 correspondence between closed polygon configurations on the left 
graph and the dimer configurations on the right graph. An edge is present in the closed 
polygon configuration of the left graph if and only if it is present in the dimer configuration 
of the right graph. Assume Fmn has weight tanh Jg on edges of T^nj and weight 1 on all 
the other edges. In other words, the edge with weight tanh Je of ^'iid the edge with 
weight e^"^^ of are dual edges. Then we have 

^T^n,-f = 2 n exp(-Je)ZT^„,Doo = n ^^^^'^^^Fmn,D 
eeE{Tmn) eeE{Tm„) 

where Zp^^ ^ is the partition function of dimer configurations on F^n- Hence we have 

ZTrr.n,Doo = J] ( ^ + eXp(2 Je) ) • 

More generally, we can expand all the Zf^„,Do^ as follows: 

ZFr.r.,De,r = U (1 + exp(2 Je) ) ((- 1)^ (- 1)^) . 

Zp ^(—1,1) is the dimer partition function of Fn with weights of edges crossed by jx 
multiplied by —1. Similarly for Zp ^(1,-1) and Zp ^(—1,-1). Therefore 

ZFmn,Doo = . max ZF^^,De,r 
6»,re{0,l} 
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Now we consider a Fisher graph Fmn- Fmn is the same graph as F^n except edge 
weights. Namely, mn has weight 1 on ah the triangle edges. Around each triangle, we have 
an even number of connecting edges satisfying 

1 

We = — , 
We 

where i(;e(resp. We) is the weight of edge e for the graph Fmn(resp. We)- All the other edge 
weights satisfy 

We = We- 

Without loss of generality, we assume that an even number of edges with weight strictly 
less than 1 are crossed by 7^;, and an even number of edges with weight greater than 1 are 
crossed by 7^. Then there is a 1-to-l correspondence between configurations in Zp^^^Dg^ 
and Zp by changing the configurations on all the edges with weight strictly less than 

1. Hence we have 

{e:We<l} 



As a result 



r) = max n (23) 

-Tmn, -1^00 Te{0 1} '^mn,L/e^T ^ > 



Proposition 4.3. Assume all the triangle edges have weight 1, and all the non-triangle 
edges have weight not equal to 1. Assume around each triangle, an even number of edges 
have weight strictly less than 1. Assume the size of the graph m and n are even, and the 
number of edges crossed by ctnd 'jy with weight strictly less than 1 are both even. Then 
P{z, —1) = have no zeros on the unit circle T. 

Proof. When both m and n are even, we have 

PfK{l,l) = Zoo - Zoi - Zio - Zu 

PfK{l,-l) = Zoo + Zoi - Zio + Zn 

PfK{-l,l) = Zoo - Zoi + Zio + Zn 

PfK{-l,-l) = Zoo + Zoi + Zio - Zu 



By (|23j), PfK{l, -1) > 0, PfK{-l, 1) > 0, PfK{-l, -1) > 0, given ah the edge weights 
are strictly positive. Hence P{z, —1) have no real roots on T. 
Assume P{z) = has a non-real zeros zo G T^. Assume 



^p = gV-i"o-^ G (0,1)U(1,2) 

If Qo is rational, namely oq = |, then after enlarging the fundamental domain to m x qn, 
Pq{z^, —1) = 0, while 2'^ is real, which is a impossible. 

Now let us consider the case of irrational ao- There exists a sequence £fc G N, such that 

lim Zq* = 1 

k—^oo 

In other words, if we assume Zq' = e^^"^"*'^ where G (—1; 1)) then 

lim Ok = 0. 
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According to the formula of enlarging the fundamental domain, 



By (23), 



P(l, -1) = (Zoo + Zoi - Zio + Ziif > (Zoi + Zii)2 



Therefore we have 



1 < lim 

fe— >oo 



(^4,01 + -^4,11)' 



(24) 



On the other hand 



lim 



^4 -^4(1,-1) 



fe.oi + ^4,11) 



lim — — 



lim 



-l7r\ 



t (2*)! 



1 61=0 



(25) 
(26) 
(27) 



Since 



0<i<m 



where P^^^*-'^ is the signed sum of loop configurations winding exactly j times around the 
cylinder (see Lemma 2.1), we have 



and 



1 



i,=o= E ^f\-^ypj''\ 

l<j<m 



lim ^ ^ 

< lim 



E 



(7rafe)2*a2*P,,(eV-i'^) 



(2t)! 



{Zi^fli + ^4,11)' 



^(vrmofc)^ V- |p(4)| 



t=l ^ ' l<j<m 

Using the same technique as described in the proof of Theorem 1.2, we have 
|pj^fe)| < Partition of configurations including nonplanar odd loops 

l<j<m 

+Partition of configurations with only even loops 
< (2^(7^"" + l)Partition of even loop configurations < C^Zj^ 

Moreover, there is a one-to-one correspondence between dimer configurations in Zoo ^-nd 
Zoi, similarly between dimer configurations in Zio and Zn. We divide the m x l^n torus 
in to i^n circles, each circle has circumference m. Fix one circle Co, and fix configurations 
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out side the circle and on the boundary of the circle. For each dimer configuration in Zqo, 
if we rotate the configuration to alternating edges along Co , we get a dimer configuration 
in Zqi. An example of such an transformation is illustrated in the Figures 11 and 12. 
As a result, 

(-^4,00 + Zif^fli - Zi^^io + Zi^^iif < C"^{Zi^fli + Zi^^iif 
where C is a constant independent of k. Hence 



lim 



fc^-OO {Zi^^Ql + Z^^ 



11; 



\2 



E 



2t fl2t p {^y^6\ 



(2t)! de^t 



\e+o 



t=i 



< hm 1 ("^^«fc)'* ^m ,72 



< 



(rmrak) 
k^i^^^-^ ^ (2t)! 



limC3™X; 



Since m is a constant, and limfc_>oo Oik = 0, we have 



which is a contradiction to (27). Therefore the proposition follows. 



Lemma 4.4. P{z,w) > 0, V(z,w) G T^. 



□ 



Proof. Assume there exists (zq, wq) £ TT^, such that P{zo,wo) < 0, since P{z, —1) > and 
P{-l,w > 0), if we consider z = e''^,w = e*"^, zq = e'^°, wq = e*"^", and {9,(t)) G [-7r,7r]^, on 
the boundary of the domain [— tt, vr] , P(e , e*"^) are strictly positive, while P(e^^o, e^'^^o) < 0, 
by continuity there exists a neighborhood O/i = ((/)o — h, cpo + h), such that for any , such 
that for any (j) S Oh, -P(e*^°, e**^) < 0. Consider straight lines in [— vr, vr]^, connecting {6q, cj)) 
to (vr, (/)), namely 

7<^(i) = {{tt -9o)t + 00,(1)) 
for G {(pQ — h,(j)Q + h), there exists t^, such that 

By continuity, there exists a (/> in the open interval {(po — h,(po + h) such that ^ a of -1. For 
instance (j) = '^■k m. reduced form where q is odd. Then after enlarging the fundamental 
domain to pra x n, we have P{ 

gi[{7r-eo)i0+9o]^ _1) = 0, which is a contradiction. 

□ 



4.4 Proof of Theorem 1.2 

In this section, we prove that the spectral curve of Fisher graphs associated to any periodic, 
ferromagnetic Ising model is a Harnack curve, and its intersection with the unit torus 
= {{z,w) : \z\ = 1, = 1} is either empty or a single real intersection of multiplicity 

2. 
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First of all, it is proved in |Dubl2j that the Ising characteristic polynomial P{z,w) 
has the same zero locus as the characteristic polynomial of the square-octagon lattice, 
denoted by Pciz, w), with suitable clock-wise odd orientation. Pc{z, w) is the characteristic 
polynomial of a bipartite graph with positive edge weights and clockwise-odd orientation, 
it is proved in jKO 06| that Pc{z,w) = is a Harnack curve, i.e, for any fixed x,?/ > 
the number of zeros of Pc{z,w) on the torus Tx,y = {{z^w) '■ \z\ = x, \w\ = y} is at most 
2 (counting multiplicities). Hence the number of zeros of Pj^{z,w) on the unit torus is 
at most 2. Since Pjr{z,w) > 0, \/{z,w) G T^, each zero Pjr{z,'w) on is at least of 
multiplicity 2. If Pjr{z,w) have a non-real zero {zq,wq) on T^, then {zo,wq) ^ {zo,wq) is 
also a zero of and each of them is at least of multiplicity 2. Then Pjr(z,w) has 

at least 4 zeros on T^, which is a contradiction to the fact that Pj^{z, w) = is a Harnack 
curve. Hence the intersection of Pj^{z, w) = with is either empty or a single real 
intersection of multiplicity 2. 

Example 4.5. (1 x n fundamental domain) Assume the Fisher graph has period 1 x n, 
and all edge weights are strictly positive. The intersection of P{z,w) = with is either 
empty or a single real node. That is, they intersect at one of the point (ibl,ibl) and the 
intersection is of multiplicity 2. 

Proof. Without loss of generality, assume P{z,w) is linear with respect to z and ^, then 
all terms in P{z,w) fall into two categories 

i) occupy z-edge exactly once 

ii) does not occupy z-edge or occupy z-edge twice 

1 X n fundamental domain is comprised of n 1 x 1 blocks, for each block, we give weights 
aii,bii,Cii,ai2,bi2,Ci2 for triangle edges, assume Oj = aiiai2, h = biibi2, Ci = CiiCi2, 1 < 
i < n. See Figure 6. Obviously such an arrangement of edge weights is gauge equivalent 
to the arrangement giving all triangle edges weight 1 and non-triangle edges positive edge 
weights. 




Figure 12: Fisher Graph on a Cylinder 

Each configuration in case i) consists of a single even loop and some double edges, the 
partition of all configurations in i is 

n 1 1 

Pi = -TT[(aj - biCi)w + (cj - aibi)] TT[(ai - ha) h (q - aibi)]z 

i=l i=l 

Configurations in case ii) depend on configurations of each 1x1 block, which are determined 
by configurations of boundary edges. Let V^^''' denote the partition at i — th block when 

2 

both edges in z direction are unoccupied, ' denote the partition at i — th block when 

2 

left edge in z direction is occupied twice, while right edge is unoccupied, similarly for V- ' 
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2 2 

and Then we have 

taO.O 2,2, / I l \ 

i^i = «j + Cj + aiCi{w + —) 

V°'^ = aiicabi2i—-w) 
w 

V^'° = biiai2Ci2{w ) 

w 

Let ki & {0,2}, 1 < i < n, then partition of ah configurations in ii is 

n 

fei ,.. . i — 1 

Assume w = e^^*^, then V^*'^ > 0, V^^'^ > 0. Periodicity implies that V^''^ and V^'^ always 
appear in pairs, then all terms in Po,2 where sin (/) appears depend only on sin^ (f>, moreover, 
given all edge weights are positive, each term with sin (f) is nonnegative, therefore 

P{Z,W) = Pl + Po,2 

n 1 " 1 

= - JJ[(ai - biCi)w + {ci - aibi)] JJ[(oj - kci) h {ci - aibi)]z 

z . w 

1=1 1=1 
n 1 " 1 

+ Yliaj + c\ + aiCi{w + -)] + JJ[6f - bi{w + -) + !]+ F{w) 

where F{w) > 0. Assume G{z,w) = z[P{z,w) — -F(w)]. Fix w, G(z,w) is a quadratic 
polynomial with respect to z. Assume 



n J 

Co = J][{ai-biCi)— + {ci-aibi)] 
i=l 

Ci = f[[a'f + c'f + aiCiiw + ^)] + fl[b'f-bi{w+^^) + l] 

i=i 

then the root for G{z, w) = will be 



w -^-^ w 

1=1 1=1 



-Ci±{Cf-4\Co\^)-2 

Zl,2 = 



2C, 

On the other hand. 



Cf - i\Cof > 4[JJ(a2 + c"^ + 2aiCi cos (^)(6f + 1 - 26^ cos 0) 

i=l 

n 

- - biCi)"^ + (cj - aA)'^ + '^{.a.i - biCi){ci - aibi) cos^]] 



1=1 



(Oj + Cj + 2aiCi cos + 1 - 26j cos 0) - [(oj - ftjCi) + (cj - 0^6,) + 2{ai - biCi){ci - aibi) cos (/>]] 

= AaibiCi s\T? (f) > 
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Therefore, Cf — 4|CoP > 0. Given ai,bi,Ci > 0, equality holds only if w is real. If 
Cf-A\Co\ > 0, then |zi,2| ^l,G{z,w) has no zeros on . Since G(l, 1) = det 1) > 0, 
P{z,w)-F{w) > onT^, so P{z,w) > onT^. If Cf-4|CoP = 0, then is real, zi = Z2 
are real, P{z,w) has real node on T^. □ 
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